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We observe that when a local patch in a radiation filled Robertson-Walker universe inflates by 
some reason, outside perturbations can enter into the infiating region. Generally, the physical 
wavelengths of these perturbations become larger than the Hubble radius as they cross into the 
inflating space and their amplitudes freeze out immediately. It turns out that the corresponding 
power spectrum is not scale invariant. Although these perturbations cannot reach out to a distance 
inner observer shielded by a de Sitter horizon, they still indicate a curious boundary effect in local 
inflationary scenarios. 



I. INTRODUCTION 



A remarkable property of the observed universe is its 
large scale homogeneity and isotropy. For instance, the 
temperature of the cosmic microwave background (CMB) 
is isotropic one part in 10^ across the sky. This smooth- 
ness appears to contradict with a generic big-bang since 
in that case one would expect to see imprints of different 
exotic objects on CMB. In particular, white holes should 
have been created in big-bang because a common big- 
crunch necessarily contains black-holes and big-bang can 
be viewed as a time reversed big-crunch. An important 
open problem in modern cosmology is to explain why the 
universe appeared in this very special state. 

It is generally claimed that the standard cosmological 
model cannot explain the isotropy of CMB, since in this 
model most of the cosmic photons we observe today orig- 
inates from causally disconnected regions in space and 
thus thermalization cannot take place to yield a uniform 
temperature. A plausible way to solve this difficulty (and 
others like flatness and monopole problems) is to modify 
the standard model by assuming an early period of accel- 
erating expansion, called inflation. Although the details 
of how the expansion is achieved depend on the model 
one considers, in general there is a scalar field whose po- 
tential energy acts like an effective cosmological constant 
yielding a de Sitter phase. By the huge exponential ex- 
pansion encountered in this period all inhomogeneities 
are smoothed out and the spatial sections are flattened 
out. The temperature drops enormously at the end of in- 
flation. However, following a reheating process the tem- 
perature again raises and one ends up with a radiation 
dominated universe. 
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Since in a cosmological model without inflation causal- 
ity precludes thermalization, it seems impossible to jus- 
tify special properties of CMB. However, as pointed out 
by Penrose (see e.g. [1]), it is paradoxical to view the 
present specialness as a feature that should be explained 
and use thermalization as the main mechanism to ac- 
count for it, since this already means that the universe 
should have been more special in the past. Namely, one 
should try to understand possible cosmic origins of the 
second law of thermodynamics instead of using it in ex- 
plaining the current state of the universe. Furthermore, 
without understanding Planck scale physics it is impos- 
sible to determine the degree of anisotropy in the ab- 
sence of inflation. Therefore it is not possible to view 
the observed isotropy as a feature contradicting with the 
standard model. 

Although inflation is proposed as a model explaining 
the homogeneity and isotropy of the observed universe, 
one usually assumes a homogeneous and isotropic dis- 
tribution of a scalar field (having a suitable potential) to 
realize it. To resolve this conflict one supposes that infla- 
tion occurred in a homogeneous and isotropic local patch 
in a larger (radiation dominated) Freedman-Robertson- 
Walker-Lemaitre (FRWL) spacetime. One then wonders 
whether inflation is 'natural' in this set up. In some cases, 
it is possible to prove a cosmic no hair theorem which 
asserts that a positive cosmological constant eventually 
dominates the cosmic evolution irrespective of initial con- 
ditions and matter content [2, 3] (see also [4]). However, 
this theorem does not say too much about feasibility of 
local inflation which is not derived by a pure cosmological 
constant. Although chaotic inflation appears to be an at- 
tractive point in simple scalar field models with suitable 
potentials [5], the issue of initial conditions for infiation 
is not yet settled (for a review see e.g. [6]). As illustrated 
in [7], without having a natural measure in the space of 
initial conditions, it is impossible to assign a probabil- 
ity for infiation. Remarkably, a measure which yields an 
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exponentially suppressed probability has been proposed 
recently in [8]. Moreover, naive arguments indicate that 
the entropy of the pre-inflationary patch is much lower 
than a patch with a typical big-bang, which shows that 
among randomly chosen initial conditions the probability 
of seeing inflation should be negligible (see e.g. [9]). On 
the other hand, there are some obstructions in embedding 
a "small" inflating region into a FRWL universe imposed 
by the propagation of null geodesies [10]. A conceivable 
way to remedy this obstruction is proposed in [11]. 

Even though the original motivations mainly arise from 
the "shortcomings" of the standard model, one of the 
main successes of the inflationary paradigm turns out to 
be its ability to produce a scale free spectrum of density 
perturbations appropriate for the formation of structure 
in the universe. According to the standard picture, den- 
sity perturbations originate from quantum fluctuations, 
which are spontaneously created in de Sitter space and 
amplified over superhorizon scales. In this paper, we 
point out that in local inflationary models there also ex- 
ists (quantum) perturbations which inevitably enter into 
the inflating region from outside, whose power spectrum 
is not scale invariant. These perturbations can spread 
out a horizon distance from the boundary and they may 
influence structure formation in that region. 

The plan of the paper is as follows. In the next sec- 
tion, we state our main argument. In section III, we ar- 
gue that due to the causal structure of local embedding 
some outside perturbations inevitably enter into the in- 
flating region and we illustrate how our main argument 
is realized in the thin-shell approximation. In section IV, 
we make a further check of scale dependence by match- 
ing the interior and the exterior modes analytically. We 
briefly conclude in section V. 



II. THE BASIC ARGUMENT 

Let us start by summarizing the well-known mecha- 
nism for the generation of scale free spectrum of density 
perturbations in inflation. Consider a scalar field (j) obey- 
ing the Klein-Gordon equation V'^cf) = in a spatially fiat 
Robertson- Walker spacetimc 



ds'^ = -df + a{tf {dx^ + dy"^ + dz^) 



(1) 



For a plane wave mode characterized by a comoving 
wavenumber k there are two important propagation lim- 
its. When the proper wavelength a/A; is much smaller 
than the Hubble radius Rr = 1/H, the mode will be- 
have like an ordinary harmonic oscillator with negligible 
damping. In that case, the ground state of the oscillator 
is a Gaussian wave function with spread given by (see 
appendix A) 



(2) 



On the other hand, when the wavelength a/k is much 
larger than Rh, the oscillator is overdamped and the 



fluctuation amplitude freezes out. 

During inflation Rh is constant but the physical wave- 
lengths increase exponentially. Therefore, a mode can 
evolve from an underdamped oscillator to an overdamped 
oscillator, but the opposite is not possible. One assumes 
that throughout inflation quantum modes were born in 
their ground states with wavelengths much smaller than 
the Hubble radius and evolve adiabatically with the 
spread given by (2) until their wavelengths become equal 
to Rh- For a mode with wavenumber k this happens 
when a/k = Rh and = H'^/{2k^). Later on, the os- 
cillator becomes overdamped, the wavelength leaves the 
horizon and the amplitude freezes out. The correspond- 
ing power spectrum P{k) = k^Acj)^. is independent of k, 
i.e. P{k) is scale free. After inflation ends, the Hubble 
radius grows more rapidly than the scale factor a, and 
thus the wavelengths created during inflation reenter the 
horizon eventually. 

Let us point out that a superhorizon perturbation is 
still propagating even though its amplitude freezes out. 
Consider, for instance, a mode in de Sitter space whose 
metric in conformal time tc < becomes 



ds' 



(Htcy^ {-dtl Jr dx^ + dy^ + dz^ 



(3) 



It is well known that the exact solution of — for 
a plane wave perturbation can be written as (see e.g. 
chapter 7 of [12]) 



H 



i{~ktc:^k.x) 



(4) 



where the overall normalization is fixed by the harmonic 
oscillator commutation relation 



(5) 



Note that (j)k = (Tk{—Htc) where ak is the canonical har- 
monic oscillator excitation. The subhorizon and super- 
horizon limits correspond to k\tc\ ^ 1 and fcjtcl ^ 1, 
respectively. Although the behavior of the amplitude dif- 
fers in the two limits, the phase dependence implies prop- 
agation along the null characteristic hypersurfaces. This 
should somehow be anticipated since these are nontrivial 
solutions of V^t/) = 0. Therefore, a mode being frozen out 
does not mean it is non-propagating, this simply signifies 
that its amplitude becomes time independent. 

In visualizing local inflation, one usually imagines that 
the inflating region is exponentially expanding into the 
ambient spacetime. Although it might be possible to pic- 
ture inflation using physical coordinates, this may create 
some confusion in subtle considerations. It is a lot safer to 
consider comoving coordinates in describing cosmological 
phenomena. In these coordinates the region itself does 
not necessarily expand. Inflation is realized by an expo- 
nentially growing metric in a (possibly fixed) comoving 
region. 

Another crucial point is that in general the inflating 
patch cannot be considered as a closed system. Although 
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FIG. 1: A spacetime sketch of local inflation in comoving co- 
ordinates. The initial conditions at time to is suitable for a 
spherical region E of radius ro to inflate in a radiation dom- 
inated FRWL space. Both regions are separated by the hy- 
persurface S. 



in some cases there may form "isolated" inflating uni- 
verses [17], we will see in the next section that even in 
this situation the inflating region is not causally discon- 
nected. Therefore, some exterior information can pene- 
trate inside. Indeed, one would expect a tendency of flow- 
ing from the surrounding region having positive pressure 
into the inflating region with negative pressure (see e.g. 
[13]). It is even possible for this flow to stop inflation, 
however we will assume that this is not the case. 

Consider, for instance, a radiation dominated FRWL 
universe. Let us assume that at time a local spherical 
region S of radius ro starts inflating for some reason (see 
figure 1), where r and t are comoving coordinates. Let 
S be the hypcrsurfacc that separates the inflating region 
from radiation dominated space. In general S can be a 
thick hypesurface, i.e. a four dimensional submanifold, 
but we assume that a thin wall approximation is valid. 
Denoting the scale factors as a/ and an, respectively, one 
can normalize coordinates such that aj(to) = afl(io) = 1 
and thus 



ai 



0-R = \ — 



1/2 



where H is the Hubble parameter in the inflating region. 
The corresponding Hubble radii are Ri — 1 /H and Rr = 
2t. Causality requires that ro should be smaller than 
the Hubble radius of the FRWL space at time to, i.e. 
2to > rg. Moreover, it is known that for stability vacuum 
energy should dominate over a region larger than the 
inflationary horizon, ro > 1/i? (seee.g. [10]). In any case 
there is no reason to expect a large hierarchy between 
these three scales, i.e. 



2tQ ~ ro 



1 



(7) 



Let us think about a (quantum) scalar perturbation 

that appeared at time t in the FRWL space whose phys- 
ical wavelength Ap is smaller than the Hubble radius 



2t > Xp. 



(8) 



FIG. 2: Left moving radial null rays of de Sitter (thick line) 
and radiation dominated FRWL space (thin line) drawn in the 
same comoving coordinates in units whore H = 1. Initially 
the Hubble constants are chosen to be equal to each other. 



This can be represented as an underdamped harmonic 
oscillator that was born in its ground state whose wave- 
function has the spread 



^4 



1 



2alk' 



(9) 



where k is the corresponding comoving wavenumber k = 
aji/Xp. The mode can be assumed to evolve adiabatically 
with (9) in the geometric optics limit, i.e. on the null 
rays of FRWL space. Moreover it cannot turn into an 
overdamped oscillator since the Hubble radius is growing 
faster than the physical wavelength. 

Assume now that (a portion of this) mode (containing 
a single pulse of one wavelength) enters into the inflating 
space at time t*. Since we are using same comoving co- 
ordinates in both regions, the comoving wavenumber k 
does not change. Therefore, the physical wavelength of 
the mode as it enters into the inflating region is amplified 

by 



A* 



(10) 



(6) Prom (7) and (8), we see that there is no hierarchy be- 
tween Xp and 1/H. Therefore, unless Xp is exponentially 



small compared to 1/H, (10) implies that 



A* > 



1 

H' 



(11) 



Consequently, this mode should freeze out in the inflating 
region. From (9), the fluctuation spectrum can be found 

as 



^4 = 



2al{t*)k' 



(12) 



which acquire a scale dependent power spectrum with 

P{k) - k^. 

Although in the above argument we assume that the 
mode enters into the infiating region at a single time t* 

and amplified immediately, this actually happens in a 
time interval proportional to the wavelength of the mode 
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(see next section). Therefore, it is important to note 
that the amphfication (10) does not occur immediately. 
Moreover, as it will be seen in the next section, it is 
not always possible to use same comoving coordinates in 
both regions. As a result, the comoving wavenumber may 
alter during crossing. As we will show, however, these 
complications do not change the main conclusion that 
outside perturbations acquire a scale dependent power 
spectrum. 

As these perturbations propagate into the inflating re- 
gion, they can spread out a comoving radial distance 1 /H 
from the boundary, which is the inflationary horizon dis- 
tance (this is due to the behavior of the radial null lines 
in de Sitter space, see e.g. figure 2). Therefore, they 
cannot reach out an observer whose comoving distance 
is larger than 1/H from the boundary. In models where 
To ^ the probability of seeing a typical observer 

(us) near the boundary is small and thus these modes 
would most likely be invisible to us to produce an obser- 
vational effect. 



III. A TOY MODEL AND THE CROSSING OF 
A SCALAR PERTURBATION ALONG THE 
BOUNDARY 

As pointed out above, since the inflating region has 
negative pressure as oppose to the exterior with posi- 
tive pressure, there should be a tendency of flowing into 
the inflating patch. This tendency can easily be seen by 
comparing radial null lines of each space pictured in the 
same comoving coordinates; it turns out that (when the 
Hubble constants are comparable) the null curves in de 
Sitter space is much more vertically aligned compared 
to the ones in a radiation dominated spacetime (see fig- 
ure 2). Therefore, information spreads a lot faster in the 
radiation region. 

Whether outside perturbations can enter into the in- 
flating space or not depends on the causal structure of 
the hypersurface S separating two regions. To determine 
S one should study the coupled dynamics of the metric 
and the scalar field driving inflation. This problem can 
be addressed in the thin-shell approximation [14] and in- 
deed the evolution of an inflating region (the false vacuum 
bubble) in a cosmological background has been studied 
well in this context (see e.g. [15-18]). In the thin-shell 
approximation, two spaces are glued over the boundary 
S describing the history of the bubble, which is a time- 
like hypersurface in both regions. The induced metric is 
taken to be continuous over S and the difference of the 
extrinsic curvatures is equal to the energy momentum 
tensor on the shell. 

Consider, for instance, a well known example from [17], 
where the evolution equations describing a spherical false 
vacuum bubble in a true vacuum region are explicitly 
solved. It is argued in [17] that at least for a range of 
initial conditions there forms "isolated" closed inflating 
universes in this setup. A solution, as seen by an out- 
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FIG. 3: The sketch of a solution for S given in [17] as seen 
by an outside observer, which yields an "isolated" inflating 
universe on the left. 

side observer, is pictured in figure 3. The outside region 
should uniquely be described by the black hole metric 
as given by the Birkhoff's theorem. Here, we see that 
the boundary is formed by a timelike curve starting from 
the past singularity at r = and extending through- 
out future null boundary. As pointed out in [17], if the 
evolution is foliated on suitable constant time hypersur- 
faces, one sees a closed inflating universe that is detach- 
ing.^ However, it is clear that the inflating region is not 
causally disconnected from the black hole region. Indeed 
all left moving perturbations created in region IV in the 
figure 3 enter into the de Sitter space and according to 
the argument presented in the previous section they ac- 
quire a scale dependent power spectrum. 

Since in the thin-shell formulation S is always taken 
to be a timclike hypersurface, it is inevitable that some 
exterior perturbations (propagating along null geodesies) 
enter into the inflating space. To explain how the main 
mechanism discussed in the previous section works in this 
context, we construct an explicit toy example in the ap- 
pendix B, where an inflating region and a flat Minkowski 
space are glued over the worldline of a spherical bubble. 
The line element in each patch is taken as 

ds^ = -dt^ + a{tf {dr^ + r^dnj) , (13) 

where a = 1 in the flat space and a = e^* in the inflating 
region. For a spherically symmetric conflguration, the 
two spheres in each region can be identified by the bub- 
ble itself. However, r and t coordinates are in general 
discontinuous. From appendix B (see (B20) and (B23)) 
we quote the equations describing the boundary: in the 
flat space S is given by 

r = — cosh(ar), t = — sinh(aT), (14) 



A critical objection to this picture is raised in [19] by indicating 
that the solution involves the past of a bifurcating Killing horizon 
which is known to be unstable. The stability issue in the thin 
shell formulation of false vacuum bubbles is also considered in 
[20]. 
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FIG. 4: The worldUne of the bubble boundary as seen by an 
observer in the flat space. 



and in the inflating region it reads 



r = 



Vq cosh{aT) 



sinh(ar) + ^ 



(15) 



Ht = ln(| sinh(ar) + 



4:Ha 



-|) -ln(aro). 



Here ro is an integration constant, cr > is the constant 
energy density on the boundary, a is given by 



a 



(7 

4' 



(16) 



and T is the proper time. To avoid any complication one 
can assume that H, a and thus a have the same order of 
magnitude. 

As pointed out in the appendix B, to solve the junction 
conditions in this simplified set up one should identify flat 
space as the interior and de Sitter space as the exterior 
patches. Therefore, these equations actually describe the 
evolution of a true vacuum bubble in an inflating space. 

The flat space trajectory (14) corresponds to the world- 
line of a particle moving with constant acceleration a 
equivalent to a hyperbola (see figure 4). To understand 
the bubble solution in de Sitter space better, we note 
that there is a singularity at r = Tg, where the denomi- 
nator of (15) vanishes. One should restrict r > since 
in that range t' > 0, where prime denotes differentiation 
with respect to r. Asymptotically as ar oo, r vq. 
The sign of r' is always negative when 4H^ — cr^ < 0. If 
4:H^ — a^ > 0, there is a turning point after which r' > 0. 
The corresponding paths are plotted in figures 5 and 6. 

Let us now determine the power spectrum of perturba- 
tions which are spontaneously created in the true vacuum 
region and passing over the inflating space. Consider a 
properly normalized right moving spherical perturbation 
of wavelength A in flat space 



1 



^ik(r — t) 



(17) 



FIG. 5: The worldline of the bubble boundary as seen by an 
observer in de Sitter space when 4_ff^ — < 0. 
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FIG. 6: The worldline of the bubble boundary as seen by an 
observer in de Sitter space when AH^ — > 0. 



where k = 27r/A. When a single pulse of one wavelength 
extending in between the null interval r — t and r — t + X 
reaches the boundary (see figure 7), it overlaps with a 
portion in the proper time interval (r, r -|- At). We as- 
sume that At -C t, i.e. the crossing time is small com- 
pared to the whole history. Therefore, the approximate 
time of passing can be defined as r* = T+Ar/2. To make 
sure that the inflating region already expanded enough, 
we focus on the perturbations entering at late times, i.e. 
e""^ > 1. Using (14), one finds that 



a 



sinh(Q!Ar/2). 



(18) 



Given any r* this equation relates At to A. 

We would like to match this pulse with a spherically 
symmetric perturbation in de Sitter space. The mode 
which has a comoving wavenumber k can be written as 



= A 



jj gife(r-tc) 



-tc + Y 
k 



where tc is the conformal time defined by 



Htr. 



-Ht 



(19) 



(20) 



Since this is not an oscillator mode spontaneously cre- 
ated in de Sitter space, there is a freedom in normaliza- 
tion which is indicated by the unknown constant A in 
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(19). Using (15), the comoving wavelength A of a pulse 
extending in between r — tc and r — tc + X can be related 
to T* and At as (see figure 8) 



Aaro 

-IT' 



sinh(aAT/2). 



Comparing now (18) and (21), we obtain 

2aVo 



k 



H 



■k. 



(21) 



(22) 



From (22), it is easy to see that if the wavelength in the 
flat space obeys 



a 



(23) 



the corresponding mode in de Sitter space satisfy 
k\tc{T*)\ < 1. Therefore, these excitations become super- 
horizon perturbations in the inflating region. Matching 
the interior and exterior scalar modes up to phase on the 
boundary one further finds 



A = 



k. 



(24) 



Due to this k dependence, (19) is not canonically nor- 
malized and the overall amplitude gives A(/)? ~ 1/fc. The 
corresponding power spectrum becomes scale dependent 
P(fc) ~ fc^, which is consistent with the results obtained 
in the previous section. 

At this point one may worry that the amplitude (24) 
is exponentially suppressed. However, this suppression 
arises from 1/r dependence of (17), i.e. this behavior 
is specific for a spherical wave created at the origin and 
does not indicate a generic feature. It is also interesting 
to note that (18) and (23) implies 



sinh(aAT/2) > tt. 



(25) 



This shows there is a minimum time width of boundary 
crossing for perturbations which become superhorizon in 
the infiating region. 

In general, one has to include a scattered left moving 

wave in flat space corresponding to reflection from the 
boundary (see next section). The reflection can be ig- 
nored if the radius of curvature of the boundary (in our 
case 1/a) is much larger than the wavelength A. From 
(23) we see that for e"^ ^ 1 there is a band of wave- 
lengths obeying both (23) and A <C l/a. Thus for these 
modes the reflection can be neglected and the above con- 
clusions must hold. 



IV. ANALYTICAL MATCHING 

In the previous section, in matching the interior and 
the exterior modes over the boundary, we have simply 

considered a pulse of one wavelength in both regions. Al- 
though this has been sufficient to illustrate our argument. 
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FIG. 7: A right moving spherical perturbation of wavelength 
A crossing the hypersurface S. 




FIG. 8: The perturbation in figure 7 exits as a spherical wave 
of comoving wavelength A in de Sitter space. 



to be more precise, one should actually glue the whole so- 
lutions. To address this problem, we first note that the 
Klein-Gordon equation can be expressed in terms of the 
derivative operator of the induced metric on S and 
the normal vector as 

(26) 

Since K is discontinues but does not involve a delta func- 
tion, this equation implies that two fields and (j)'^ , 
which are obeying V^t/)^ = in respective regions, can 
be glued on the boundary if they satisfy 



(27) 
(28) 



We now use these junction conditions to study the cross- 
ing of a perturbation. 

In our case, there is an incoming wave which is the 
properly normalized, right moving, spherical perturba- 
tion 



1 gife(r-t) 



2k 



(29) 



The reflected and the transmitted modes can be taken as 
the general, spherically symmetric, left and right moving 
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waves in flat and de Sitter spaces, respectively. These 

can be written as 

g(l/(r + t)) f'{r-Qt, , /(r-y 

0re/ = , 9tr = 1 , 



(30) 

where /, g are arbitrary functions and prime denotes dif- 
ferentiation with respect to the argument. The general 
solution in de Sitter space can be obtained from (19) by 
superposition. Two conditions (27) and (28) can be used 
to determine two unknown functions / and g. Note that 

4>~ = i'in + 4>reS and (/)+ = (j)tr ■ 

We consider the late time matching of the solutions, i.e. 

we assume e"'^ ^ 1. From (14) and (15), the components 
of the normal vector can be found in flat space as 



n = n 



2 ' 



while in de Sitter space they read 



AHa 



n = 



2roa^ 

IT 



(31) 



(32) 



where we only keep the leading order contribution in an 
expansion in e~"^. We define a (dimensionfuU) variable 
X = ae~°''^. In units where a = 1 we have a; <C 1. To 
first order in x, the continuity of the modes (27) implies 



xe 



ikx/a 



+ 2xg=-f--xf, 



ro 



H 



(33) 



where prime denotes differentiation with respect to the 
argument and 

g = 9{x), f = /(ro + ^x). (34) 

On the other hand, the continuity of normal derivatives 
(28) gives 




X — 




kxj cx^ 



2axg' — 2ag 



(35) 



Solving g from (33) and using it in (35) we obtain the 
following second order equation for /: 



2ro(CT + 2Q) _ 2q 2^ 
^2 ^/ ^ roH^ 



2a 



2k 



2i 



2k 



a 



-X e 



ikx I 



.(36) 



The right hand side of this equation is the source term 
related to the incoming wave. 

It is possible to solve (36) exactly. The homogenous 
solution involves two Bessel functions and the particular 

solution can be obtained by variation of parameters. We 
need to keep the unique particular solution supported by 



the source incoming wave. However, (36) has corrections 
in powers of x, i.e. both the fimctions multiplying / and 
the source in the right hand side are just the leading 
order contributions. Therefore, only the first term in the 
expansion of the exact solution in x can be trusted. 

Since the source term in (36) is oscillating, the solution 
can be written as / — e*'^^/" Fix), i.e. / will oscillate 
exactly the same way as the source. From (34) we then 
find 



f{y) 



JkHy/{2a^ro) 



(37) 



thus the comoving wavenumber k in de Sitter space can 
be identified as A: = kH/{2a'^ro), which is precisely the 
relation (22) derived in the previous section. 

Let us now determine the solution of (36) in the short 
and long wavelenght limits. To avoid any complication 
we assume that H, a (and thus a) and 1/ro have the 
same order of magnitude. 

We start with the short wavelength limit kxjo? ^ 1, 
where the first term in the right hand side of (36) can be 
neglected and the solution can be expanded as 



ikxjo? 



(l + 0(x)). 



(38) 



In (36), /" term becomes much larger than the two other 
terms in the left hand side and the constant A can be 
fixed as 



(o- + 2a) fcV2fc' 
Here k dependence is crucial. Using / in (30), we find 



(39) 



^^^-^Hk|2c? 

(cr + 2a) 



U gik{r-t^) 



2k 



-ic+~). (40) 



The constant in the square brackets is an order unity con- 
stant. Comparing with (4) and nothing that k and k are 
nearly equal to each other (recall that we assume H, a 
and 1/ro have the same order of magnitude) , we see that 
the transmitted wave is a properly normalized mode in 
de Sitter space. Indeed, in the limit that we consider 1/k 
term can be neglected compared to tc in the parentheses, 
which shows that this is a subhorizon mode. Therefore, 
perturbations which have short wavelengths in flat space 
can be matched by subhorizon modes in the inflating re- 
gion which acquire a scale free power spectrum. This is 
consistent with the results of the section II since the am- 
plification (10) is not enough to push the physical wave- 
lengths of these modes out of the inflationary horizon. 

Consider now the long wavelength limit kxjo? <C 1. 
In this limit the second term in the right hand side of 
(36) can be ignored, the exponential can be set to unity 
and the solution can be expanded as 



j = A{\ + 0{x)). 



(41) 



The amplitude can be flxed from (36) as A = rf)H/\/2k 
which yields a spherical, superhorizon mode in de Sitter 
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space 



2kr 



(42) 



However the normalization is different than the normal- 
ization of the canonical superhorizon perturbation, which 
is l/k^/2k . This gives a scale dependent power spectrum 
with P{k) ~ k'^, which is consistent with the results pre- 
sented in the previous sections. 



V. CONCLUSIONS 

Inflation offers a successful way of generating a scale 
free spectrum of density perturbations consistent with 
cosmological observations. According to the general 
view, all classical inhomogeneities are smoothed out by 
the exponential expansion, leaving only room for quan- 
tum fluctuations. Naturally, quantum modes are as- 
sumed to be born as harmonic oscillators in their ground 
states, which have physical wavelengths smaller than the 
inflationary horizon. The wavelengths are then pushed 
to superhorizon sizes by expansion and their amplitudes 
freeze out. After inflation ends, they reenter the horizon 
in a radiation dominated phase as classical perturbations 
seeding cosmic structure. 

The success of the above mechanism profoundly de- 
pends on the very special properties of the generation 
and the propagation of quantum fluctuations in de Sit- 
ter space. Indeed, even in de Sitter space the well known 
ambiguity in the choice of an invariant vacuum may alter 
some of the predictions of the theory (sec e.g. [22]). In 
this paper, we point out that in local inflationary models 
the causal embedding of an inflating patch into an ambi- 
ent spacctime generally allows outside pertiirbations to 
enter into the inflating region. It is clear that these per- 
turbations may not have a scale-free power spectrum. 

In this paper we show that quantum fluctuations en- 
tering into a local inflating patch have a scale dependent 
power spectrum. In section II, we first give a general 
argument and in sections III and IV we illustrate it by 
an explicit toy example in the thin-shell approximation. 
Although the set up studied in this context is artificial 
and cannot be considered as a part of a realistic scenario, 
it nicely touches to the salient points raised in section II. 
It would be interesting to construct a numerical or an 
asymptotic false vacuum bubble solution in a radiation 
dominated FRWL space and verify that scale freeness is 
spoiled for the quantum modes crossing into the inflating 
region, as suggested by the basic argument of section II. 

One may think that the above conclusions are avoided 
in models where the inflating patch detaches from the 
ambient spacetime. We consider one such example from 
[17] and indicate that the patch is still causally connected 
to the surrounding region. Indeed, in the thin-shell ap- 
proximation causal detachment is impossible since the 
boundary in between the two regions is taken to be a 
timelike hypersurface in both sides. 



The key issue is then to determine whether causal dis- 
connection is possible in local inflationary models at all. 
One of the nicest features of infiationary paradigm is that 
it uses classical general relativity in determining the evo- 
lution of the metric and thus its predictions are not sensi- 
tive to the quantum gravitational corrections. In general 
relativity, on the other hand, causal detachment implies 
the existence of a real event horizon. To our knowledge 
black holes are only examples where observer indepen- 
dent causal separations take place. In some models using 
chaotic inflation, it is argued that after big bang (or after 
universe emerges from a spacetime foam) regions suitable 
for infiation expands and others rccoUapse. In such local 
models it seems that the above conclusions are avoided. 
However, it is not clear to us how this picture can be 
verified without understanding the end point of a col- 
lapse, which requires a detailed knowledge of quantum 
gravity. Indeed, even in semiclassical gravity black holes 
radiate energy and thus it appears that a complete causal 
detachment is not possible. 

Finally, these flndings do not imply a disagreement 
with observations in local inflationary models. To claim 
a discrepancy, one should first make sure that the con- 
tribution of outside quantum perturbations to structure 
formation is not negligible. Moreover, as pointed out 
previously, depending on the size of the initial infiating 
patch and the location of the observer inside, the out- 
side perturbations may not be able to reach the observer 
since in de Sitter space there are particle horizons and an 
observer is shielded by a horizon. In any case, our results 
indicate that the boundary effects should be considered 
as an issue in local inflation. 



APPENDIX A: SCALAR FIELD QUANTIZATION 
IN A FRWL SPACE 

In this appendix we review the quantization of a real, 
massless scalar field in a FRWL spacetime with the met- 
ric 



ds'^ = -dt^ + a{tf {dx^ + dy"^ + dz^) 
The canonical action can be taken as 



-gg^''d.cj)d^(l), 



(Al) 



(A2) 



which gives V^^ = 0. Since in a free field theory different 
Fourier modes decouple, we consider a single excitation 
with a fixed comoving wave- vector k: 



where the field equations imply 



4>k + 3H<pk 



T<t>k = 0. 



(A3) 



(A4) 
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From (Al) and (A2) the momentum conjugate to (f> be- 
a^(p. The quantization can be achieved by 



comes 
imposing 



|0>=0, 



(A5) 



where |0 > is the vacuum as seen by a comoving observer. 

For a subhorizon mode k/a^ H the friction term in 
(A4) can be ignored and an approximate solution for 
can be found as 



t)k = aue 



ikt/a 



The commutator in (A5) then impUes 



1 



2ko?' 



(A6) 



(A7) 



Therefore, ak and aj^ arc creation and annihilation oper- 
ators and the system is equivalent to the harmonic oscil- 
lator. One can easily calculate the two-point function or 
the spread of the Gaussian ground state wave-function 

as 



By Einstein's field equations the discontinuity of the 
extrinsic curvature can be related to the energy momen- 
tum tensor on the hypersurface Sij. Moreover, the four 
dimensional energy-momentum conservation also implies 
a conservation equation for Sij . One can show that these 
respectively yield the following equations for the bubble 
(known as junction conditions), 



[-Kij + Khij\t = Si 



+ 



(B3) 
(B4) 



where Di is the covariant derivative of hij. 

In cither region, the tangent vector corresponding 
to the derivative operator dr can be calculated as 



: t'dt + r'dr, 



(B5) 



where prime denotes differentiation with respect to the 
proper time t. By definition v^^Vfj, = — 1. The two tan- 
gent vectors on the sphere together with w** form a bases 
in the tangent space of S. This information is enough to 
determine the unit normal vector as 



A<i>i =< o|<A^|o >= 



(A8) 



APPENDIX B: DERIVATION OF THE BUBBLE 
SOLUTION 

The aim of this appendix is to obtain an explicit equa- 
tion for the hypersurface S which is separating a fiat 
Minkowski space from an inflating region. We start by 
considering the actual problem, i.e. a false vacuum bub- 
ble in a FRWL space. The metrics in both regions can 
be written as 



-dt\ 



-a±{t±f {dr 



[dni 



(Bl) 



where -|- and — denote the exterior and interior patches, 
respectively. For a spherical bubble one can identify the 
interior and the exterior spheres with the bubble itself. 
However, t and r coordinates are in general discontinuous 
over 5* and they should be differentiated in both regions. 
The third coordinate on S can be chosen as the proper 
time r. The hypersurface will be completely specified 
when the trajectories t±{T) and r±{T) are solved. By 
spherical symmetry, the induced metric hij on S, where 
i,j indices refer to the coordinates (r, 6, (j)), can be writ- 
ten as 



(B2) 



In the thin-shell approximation hij is assumed to be con- 
tinuous over S. Note that R{t) is the proper radius of 
the bubble at time r. If one is only interested in the in- 
trinsic properties, it would be enough to determine R{t). 
However, wc arc also concerned with the embedding of 
the bubble into the ambient spacetime. 



= {ar'dt+a~^t'dr) sign(t') 



(B6) 



where the sign function is inserted to make sure that is 
pointing from inside out, i.e. in the growing r direction. 
Since S is now uniquely specified, it is straightforward to 
calculate the extrinsic curvature which has components 



dt 



dt 



Kee=(a'r''^r' + art']sign{t') 
K^^ = Keel siv? 9, 



sign(^'), 
(B7) 



where the last relation follows from spherical symmetry. 

For a spherical false vacuum bubble associated with 
a minimally coupled scalar field, the energy momentum 
tensor on the hypersurface can be assumed to have the 
form [17] 



Sij — dhij . 



(B8) 



Thus in our problem there are five dynamical fields of 
interest which are t±{T), r±{T) and cr(r). The fact that 
T is the proper time gives two equations 



t'l - alr'l = 1. 
The continuity of the induced metric yields 

d2 2 2 2 2 

R = a+rz^ = a_r_, 



(B9) 



(BIO) 



where R is defined in (B2). The conservation of stress 
energy tensor (B4) can be satisfied provided 



a' = [at'r' (p + P)] + 



(Bll) 
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Finally {09) component of (B3) implies 



(a^ r'^^r' + art') 



--a r 



(B12) 



where we assume sign(t') > in both regions. Equations 
(B9), (BIO), (Bll) and (B12) supply enough information 

to determine the dynamics of the bubble. Let us point 
out that (rr) component of (B3) gives a second order 
equation 



ar't" - at'r" + , 



> da 
■—1 
dt 



da 
dt 



a 

r 



(B13) 



which is satisfied identically provided the above five equa- 
tions and the Einstein's equations in the exterior and 
interior regions hold. 

Following [21] it is possible to obtain a useful equation 
for R. After some simple algebraic manipulations one can 
show that {99) component of the discontinuity equation 
(B3) implies 

Kj-Kti + ^\ . (B14) 



From (B7) it follows that K^g can be expressed in terms 
of R as 



K? 



06 



i?2 [l+i?'2_jy2^2] 

Using (B15) in (B14) yields 



(B15) 



R^ 



2„2 



R'a 



R" + i = ^{Hl-Hl) +-{Hl + Hl)+ , 

(B16) 

where H± = dlna±/dt± are the corresponding Hubble 
parameters. Although this looks like a single equation for 
R, let us note that there is a dependence on the variable 
a and a hidden dependence on through . Note that 
H- = H is constant in the interior de Sitter space. 

It turns out that it is difficult to solve the above field 
equations to obtain an explicit solution for the boundary. 
Since our aim in this paper is to examine how perturba- 
tions propagate from one region into the other, rather 
than studying the bubble dynamics, we make some sim- 
plified assumptions. 

First note that for local inflation it is desirable to have 
iJ_ ^ Therefore, in (B16) one can ignore the terms 
involving . This is nearly equivalent to taking outside 
region to be the flat space, except a is still time depen- 
dent. In (Bll), the interior contribution to the right 
hand side vanishes since the cosmological constant obeys 
P- + P- = 0. By field equations, the exterior contribu- 
tion p+ + P+ is proportional to dH+/dt+ = —2/t\. We 
assume that inflation occurs at a sufficiently later time 
after big bang (i.e. t+ 1) such that this term is very 
small and thus a' ^ 0. This effectively means that the 
outside region can be taken as the flat space. Summariz- 
ing we take 

a-=e"*-, a+ = 1, ct' = 0. (B17) 
These assumptions may not be suitable for a realistic 
scenario using inflating bubbles in a cosmological back- 
ground, but they are not harmful for our purposes. 



After these simplifications, (B16) can be solved for R 



as 



where 



R = — cosh(aT) 
a 



a 



a 4' 



(B18) 



(B19) 



Since a+ = 1, (BIO) and (B9) can be used to determine 
r+ and t+ as 

r+ = — cosh(QT), f+ = — sinh(QT). (B20) 
a a 

To solve the interior fields we first note that (BIO) can 
be used to fix t- in terms of r_ as 



Ht_ =lni?-lnr_. 



(B21) 



Using (B21) in the proper time equation (B9), one gets 
a quadratic equation for In r'_ which can be solved alge- 
braically as 



_ lnfi'±gfi|^-f| 
- 1-mR^ 



(B22) 



Although it looks complicated, this equation can be inte- 
grated to obtain r_(T). Further using (B21) the interior 
embedding coordinates can be fixed as 



r_ = 



ro cosh(aT) 



AHcr 



I sinh(ar) ± 

4H^ — 

= ln(| sinh(aT) ± — — 1) - ln(Q;ro). 

4/7 c 



(B23) 



The correlation between ± signs in (B23) and (B22) de- 
pends on the sign of 4iJ ^ — ct^ . It may first be seen sur- 
prising that there appears two solutions since the initial 
system of differential equations are first order. However, 
it is easy to see that ± solutions are related by time re- 
versal r ^ — T. Without loss of any generality we choose 
+ sign in (15). 

At this point we recall that although (B16) was ob- 
tained from the original jimction condition, the two equa- 
tions are not equivalent. Therefore, a final check is re- 
quired to see whether (B20) and (B23) obey (B12). It 
turns out that to satisfy (B12) one should either choose 
cr < or interchange the interior and the exterior re- 
gions. Instead of choosing tr < 0, which amounts to 
assume a negative cosmological constant on the shell, we 
simply switch the roles played by the interior and exterior 
spaces. 
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